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Abstract 

This paper is the first part of a three part tutorial on multidimensional generalized 
functions (GFs) and their applications in aeroacoustics and fluid mechanics. The 
subject is highly fascinating and essential in many areas of science and, in particu- 
lar, wave propagation problems. In this tutorial, we strive to present rigorously and 
clearly the basic concepts and the tools that are needed to use GFs in applications 
effectively and with ease. We give many examples to help the readers in understand- 
ing the mathematical ideas presented here. The first part of the tutorial is on the 
basic concepts of GFs. Here we define GFs, their properties and some common 
operations on them. We define the important concept of generalized differentiation 
and then give some interesting elementary and advanced examples on Green's func- 
tions and wave propagation problems. Here, the analytic power of GFs in applica- 
tions is demonstrated with ease and elegance. Part 2 of this tutorial is on the diverse 
applications of generalized derivatives (GDs). Part 3 is on generalized Fourier trans- 
formations and some more advanced topics. One goal of writing this tutorial is to 
convince readers that, because of their powerful operational properties, GFs are 
absolutely essential and useful in engineering and physics, particularly in aeroacous- 
tics and fluid mechanics. 
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Nomenclature 

C°° : The space of infinitely differentiable functions 

C n : The space of continuously differentiable functions up to order n 

D : The space of infinitely differentiable functions with compact (closed and 
bounded) support 

D' : The space of continuous linear functionals on the space D, i.e., the space of 
Schwartz GFs 

GD : Generalized derivative 
GF : Generalized function 

OF : Ordinary function that is locally Lebesgue integrable 
R : The unbounded real line 

R 2 , R 3 : The unbounded two and three dimensional spaces, respectively 
e : Belongs to 

[0, 1] \ [a, b] : Set subtraction, i.e., [0, a) (J (b, 1] 

1- Introduction 

This paper is the first part of a three part tutorial on multidimensional generalized 
functions (GFs) and their applications in aeroacoustics and fluid mechanics. As the 
name implies, GFs are the generalization of the ordinary functions (OFs) that are 
defined as locally Lebesgue integrable functions. Here, mathematical objects that 
are not OFs, such as the well known Dirac delta "function", are added to OFs and 
the resulting space of objects is called GFs. In this tutorial, we will present rigor- 
ously and clearly the basic concepts and the tools that are needed to use GFs in 
applications effectively and with ease. We give many examples from aeroacoustics 
and fluid mechanics to help readers in understanding the mathematical ideas pre- 
sented here. 
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The first part of the tutorial is on the basic concepts of GFs. Here we define GFs, 
their properties and some common operations on them. We define the important 
concept of generalized differentiation and then give some interesting elementary 
and advanced examples on Green's functions and wave propagation problems. 
Here, the analytic power of GFs in applications is demonstrated with ease and ele- 
gance. Part 2 of this tutorial is on the diverse applications of generalized derivatives 
(GDs). Part 3 is on generalized Fourier transformations and some more advanced 
topics. One of the goals of writing this tutorial is to convince readers, that because 
of their powerful operational properties, GFs are absolutely essential and useful in 
engineering and physics, particularly in aeroacoustics and fluid mechanics. 

The linear wave equation and its counterpart in the frequency domain, the 
Helmholtz equation, are the most important partial differential equations of aeroa- 
coustics. For example, both the Lighthill 1 jet noise equation and the Ffowcs 
Williams-Hawkings (FW-H) equation 2 are linear wave equations. We will primarily 
be concerned with the wave equation here. Although this equation in one, two, 
three and higher space dimensions has been analyzed since the eighteenth century 3 , 
there are still many issues and subtleties associated with finding its solution, particu- 
larly in aeroacoustics, that are in need of further attention. It is true to say, and this 
fact is well-known to mathematicians, that the appropriate mathematical tool for 
solving PDEs is the theory of distributions or generalized functions (GFs) 3-12 . The 
derivation and the solution of the FW-H equation depend heavily on the concept of 
GFs, which requires much mathematical maturity from users. Furthermore, even 
when we use the Lighthill equation, we may be dealing with flowfields that have 
discontinuities within them, such as shocks, which can be sources of sound. The 
natural tool to work with these discontinuities is GF theory. The principal goal of 
this tutorial is to give a working knowledge of multidimensional GF theory to the 
readers, particularly when applied to the linear wave equation. 

There are many books and papers on the mathematics of the wave equation and the 
related problems that are used by acousticians. Unfortunately, almost all of the 
commonly referenced books and papers were written before, or without the use of, 
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the theory of distributions (generalized functions) of Laurent Schwartz 4 . These 
books and papers use an intuitive approach to Dirac delta function and its algebra 
that is not rigorous and is prone to giving erroneous results. Furthermore, much of 
the true power of GF theory is not apparent from the intuitive approach. The avail- 
able books on mathematics that discuss generalized functions (GFs) rigorously 4- 6 
are also beyond the reach of most researchers because of the level of abstraction 
involved to comprehend the theory. We believe that most of the abstractions are not 
necessary for application of the GFs. The situation is much like the use of the real 
numbers by the man in the street in everyday business without ever knowing about 
Cauchy convergent sequences because he only needs to know the useful properties 
of the real numbers. For mathematically minded persons the abstractions in the 
development of the theory can be learned most easily after a working knowledge of 
GFs is attained. 

We will start by giving the definition of GFs as continuous linear functionals on a 
given space of functions. We will do this in the simplest possible way. We will then 
discuss how the definition of many useful operations on ordinary functions can be 
extended to the GFs. Of particular interest to us here is the definition of generalized 
derivative and its many useful applications. These include Green’s functions of 
ODEs and PDEs, finite part of divergent integrals, the Kirchhoff formula for mov- 
ing surfaces and deriving the jump conditions from conservation laws. We will 
discuss some questions associated with the derivation of the FW-H equation. In 
particular, we answer the question of whether any function other than the null func- 
tion can be assumed inside the moving surface in the derivation of the FW-H equa- 
tion. We will see that in every aspect of working with the linear wave equation, the 
natural mathematical tool is the GF theory. 

With this tutorial article, we wish to convey to readers our enthusiasm for working 
with GFs as the most essential, useful and powerful tool for application in acoustics, 
fluid mechanics and other areas of physics and engineering sciences. 

2- Generalized Functions (GFs) 

In this section we will start by giving the definition of GFs. We then present the 
most common operations on GFs. We will pay a special attention to generalized 
differentiation, which is the fundamental concept for obtaining the Green’s func- 
tion, finding the finite part of divergent integrals and working with jump conditions 
from conservation laws. We will then give many applications of GFs relevant to 
acoustics. We will say little about generalized Fourier transformation because that 
subject will be covered in Part 3 of this tutorial. 
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In learning GF theory, you should know that there are several different approaches 
to the definition of GFs. The method we present here is the functional approach of 
Schwartz 4-11 . The sequential approach is used by Temple 13 and Lighthill 14 , Jones 12 

and Mikusinski 15 . The completeness theorem of GFs relates the two approaches 7 . 
The approach by Bremermann defines GFs as the boundary value of analytic func- 
tions 16 .This subject is valuable in studying generalized Fourier analysis and, in its 
full generality, requires the knowledge of the highly specialized topic of functions 
of several complex variables. Most of Bremermann’ s book 16 does not require this 
topic and for learning about one dimensional GFs, it is easily accessible to readers 
with a good knowledge of analytic functions of a single complex variable. Eventu- 
ally, when you work on advanced Fourier transform theory, you will have to learn 
Bremermann’ s approach because it sheds light on the derivation of many of the 
results obtained by physicists and applied mathematicians using ad hoc and mostly 
intuitive reasoning. The subject can be shown to be also closely related to the func- 
tional approach we use here but we will not say any more about this connection 7 . 

The Polish mathematician Jan Mikusinski introduced one more method of defining 
GFs known as the operational approach which is based on the concepts from 
abstract algebra 17 ’ 18 . Mikusinski’s publications are beautifully written and are easy 
to read and comprehend. Apart from explaining rigorously the Heaviside opera- 
tional calculus, this approach can also be used to solve an entirely different set of 
problems than the first three approaches, such as finding the exact solutions to finite 
difference equations and recursion relations. The relationship between Mikusinski’s 
operational approach and the functional approach cannot be established easily and 
we are not aware of any publications on the subject. We encourage readers to study 
this approach because of its usefulness and power in solving some important prob- 
lems of physics and engineering 17-20 . We will say no more on this approach here. 

And finally, readers should be aware of the fact that in recent years GF theory has 
been extended to the so-called nonlinear GF (NLGF) theory where some of the 
limitations of Schwartz distribution (GF) theory are removed 21-23 . For example, it 
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can be shown that one cannot extend the operation of multiplication to all of 
Schwartz GFs. This shortcoming of the Schwartz theory does cause some difficul- 
ties in applications that Colombeau 21 , Rosinger 22 and Oberguggenberger 23 as well 
as other mathematicians have tried to resolve. Currently, the NLGF theory is very 
abstract and difficult to master. The subject is under intense development but is 
badly in need of simplification. However, it appears that the subject of nonstandard 
analysis (NSA) 24 will eventually be able to clarify and simplify NLGF theory consid- 
erably and bring the theory within the reach of physicists and engineers 25 . We will 
not cover this subject here. Readers should know that to understand NLGF theory 
and NS A, requires one to have a knowledge of several esoteric areas of mathemat- 
ics including topology, abstract algebra and mathematical logic. Fortunately, there 
are many expository books and articles on these subjects on the internet, particu- 
larly on NS A. We will not cover this theory here. 

For those interested in the history of GF theory, we recommend the books by 
Liitzen 26 and Dieudonne 27 as well as the article by Synowiec 28 . In the opinion of 
Liitzen, the Russian mathematician Sobolev invented the concept of GFs (more 
exactly, generalized derivatives) and Schwartz developed the theory of GFs without 
being aware of Sobolev’s earlier contributions 4 . As expected, many other mathemati- 
cians contributed to the field. Without a doubt, the members of the Moscow School 
of Mathematics under I. M. Gelfand have contributed substantially to the develop- 
ment and applications of GF theory making the subject accessible to physicists and 
engineers. The well-known book on GFs by Gelfand and Shilov 7 is the first of a six 
volume opus by this School published in the West in the sixties of the last century. 
This monumental work, which is still read today, has popularized the GF theory 
enormously among physicists and mathematicians by its breadth and depth as well 
as by the clarity of the exposition. For learning about more recent advances we 

suggest the books by Hormander 5 and Taylor 6 which are both written for 
mathematicians. 
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Most of this section is based on the first author’s expository NASA publication 29 . In 
addition to this publication, we recommend that you consult the books by 
Strichartz 10 , Stakgold 8 , Gelfand and Shilov 7 , and Kanwal 11 . For beginners, we partic- 
ularly recommend the excellent book by Strichartz for deeper understanding of 
many of the concepts discussed here. This book emphasizes GFs in one dimension 
which is essential in some applications, e.g., time series analysis. However, there 
are many problems of physics and engineering for which multidimensional GF 
theory is very important 7-9 ’ 11,12,29 This is the subject that we concentrate in the 
present paper. The emphasis will be on using the theory to solve problems rather 
than stating theorems and giving proofs. 

■ 2.1- Definition of Generalized Functions (GFs) 

What we call a function, say fix) , in calculus can be thought of as a table of ordered 
pairs (x, fix)). We can visualize this table as a graph of x versus y - fix). If we 
want to generalize the concept of a function, we must change the definition of a 
function in such a way that our old functions are still considered functions and we 
have many new useful mathematical objects that are functions by the new defini- 
tion. Once we introduce a new definition of a function, we must make sure that we 
can extend all the useful operational properties of the old functions, e.g., addition of 
two functions, to all the new mathematical objects which we now collectively (the 
old functions and the new objects) call generalized functions (GFs). For now, we 
refer to all the functions that we learned about in calculus as ordinary functions 
(OFs). We will elaborate further on this below. 

Let us start by considering an ordinary function g(x) which is periodic over the 
interval [0, 2 tt ]. We know that the complex Fourier components of this function are 
given by the relation 

1 r 2n 

G n - — I g{x)e inx dx (/? = — ootooo) 

In Jo 


( 1 ) 
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We know that g(x) = Yj7i=-oo G n e inx , that is, the periodic function g(x) can be con- 
structed from the knowledge of the table with infinite entry {G n ,n = -ootooo}. 
From another point of view, the set of complex numbers {G n } can be considered as 
the mapping of the space of functions [e~ in x , n - -oo to ooj into complex numbers 

by the rule given in eq. (1). Such a mapping from a given space of functions into 
scalars (real or complex numbers) is called a functional. Remember that a func- 
tional is always a rule that describes the mapping from a given space of functions 
into scalars. The numbers {G„} are the functional values of this functional. Note 
that the function g(x) enters the definition of the functional. As a matter of fact, we 
should use the functional notation G [e~ 1 n x ] = G n to remind us that the numbers 
G n , n — —oo to oo depend on the function g(x). So far we have shown that the peri- 
odic function g(x) can be described (identified, defined) by the table of its func- 
tional values: 


(G n = G[e~ inx ] = ^- f 2n g(x)e- inx 
K 2 71 Jo 

on the space [e~ inx , n - -ootooo} 


n — — oo to oo 


( 2 ) 


We have assumed here that the function g(x) has a Fourier series (e.g., satisfies the 
Dirichlet conditions, etc.). Now we have the germ of an idea to extend the defini- 
tion of a function. 


Let us take a space of functions of one variable D with a large number of well- 
defined functions in the space. We will say more about the properties of the func- 
tions in this space soon. Let us select an ordinary function f(x) and for any function 
0 (x) in the space D, define the functional (i.e., the rule): 

/'XX) 

F[(j)\ = I f{x) 0 (x) d x (/ is fixed in this rule and 0 varies all over D ) (3) 

kJ — OO 

Note carefully how the function /(x) enters the definition of the functional. To 
describe /, this function is kept fixed in the mle given by eq. (3) and different func- 
tions 0 are taken {conceptually , of course, not literally ) from the space D to form 
the table {F[<fi], 0 eD}. Therefore, we must have so many functions in space D that 
any two different ordinary functions generate different tables of functional values at 
least for some entries. The space D is known as the test function space. The ques- 
tion is how to select the space D in such a way that the set of the functional values 
{F[<p], 0 e D\ represents (i.e., describes or has as much information content as) any 
ordinary function /(x). Before we do this, we must now elaborate on what we mean 
by ordinary functions. 
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By ordinary functions (OFs), we mean the locally Lebesgue integr able functions. 
This means that they are functions that are Lebesgue integrable 7- 11,19 over any 
finite interval (see Note 1 for some comments on Lebesgue integration). These 
functions include any conceivable functions, no matter how wild, that are needed in 
applications. If you do not know Lebesgue integration, you can think of Riemann 
integration that you learned in calculus courses. Lebesgue integration is a generaliza- 
tion of the Riemann integration that can be used for the integration of most func- 
tions that do not have a Riemann integral. 

One can show that there is an uncountable number of OFs. Therefore, the space of 
test functions must contain an uncountable number of functions so that any two 
different OFs can be distinguished from each other when viewed as tables of func- 
tional values. This means that our tables of functional values used to represent OFs 
are also uncountable. Furthermore, since OFs are as wild as they can be, we require 
that the functions in space D to be as well-behaved as they can be so that the inte- 
gral in eq. (3) is always defined. This leads us to require that the functions in D to 
be infinitely differentiable (denoted C°° functions). Also to ensure the integral in eq. 
(3) does not become infinite, we assume that all functions in space D are identically 
zero outside of a finite interval. This interval can be different for different test func- 
tions . A function that is identically zero outside some finite interval is said to have 
bounded support. Thus, we have come to the following conclusion regarding the 
space of test functions D that we need to use to identify OFs by their tables of func- 
tional values: 

The test function space D that one needs to use in 
defining OFs by their tables of functional values is the space of 
all C°° functions with compact ( closed and bounded ) support. 

Let us say something about this space next. The best known, and in many ways the 
fundamental, example of a function in D is the following bell-shaped function with 
the peak value at the origin. It is defined for any real number a > 0: 




\x\ ^ a 


( 4 ) 


10 
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Figure 1- The Graph of the function 0 (x; a) 

It is easy to see that the support of this function is [- a, a] and it is infinitely differ- 
entiable on the open interval (- a , a). Figure 1 shows the graph of this function. 
The proof of infinite differentiability at points x - - a and x - a is much harder 
and is an exercise in calculus albeit requiring much technical skill. You should 
know that at these two points all the derivatives of this function are zero. Are there 
any other types of functions in space D besides 0 (x; a)l The answer is affirmative 
and we will turn to this point next. 

One proves rigorously in real analysis that there is an uncountably infinite number 
of continuous functions. For any continuous function g(x) defined on the bounded 
interval I = [b, c], the convolution of g(x) and 0 (x; a) on I gives a C°° function 
with the bounded support [b - a, c + a ] 19 '. Therefore, there is an uncountable num- 
ber of functions in the test function space D as we are required to have. Another 
fact to remember about this space is that the functions in space D are not analytic 
even though they are C°° functions ! 

At this stage, we have a new way of looking at OFs but our original intention was 
to extend the definition of a function so that we have new objects that are not OFs 
but are considered functions. To do this we study two important properties of the 
functionals described by eq. (3). The first one is very easy to see: 


1 - The functional F[0] is linear. 
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This simply means that for a and ft scalars and 0 and if/ functions in D, we have 

F[a 0 + (3 f] - a F[f\ + /3 F[if/] (5) 

To show this, replace 0 (x) on the right of eq. (3) with a 0 + fif and use the linear- 
ity property of integrals followed by functional interpretation of the result to get the 
right side of eq. (4). We remind the readers that the OF /(x) used in the definition 
of the functional of eq. (3) is fixed here and the functional F[f] is, in point of fact, 
representing fix) . 

The second property is much harder to explain and comprehend. It is: 

2 - The functional F[f] is continuous. 

The first thing you must know is that this property, whatever it means, gives GFs a 
lot of useful properties that we need in applications. Second, to beginners the defini- 
tion of continuity of a functional appears too complicated and perhaps a futile 
attempt in abstraction. We admit that the definition of the continuity of a functional 
is complicated but its definition cannot be stated in simpler form than we give here. 
The original definition by Schwartz was much more complicated and based on 
locally convex topological spaces 4 . Just remember that you will appreciate the full 
usefulness of this property when you learn some more about GFs and their applica- 
tions. See, for instance, Example 3 below. 

We first distinguish some sequences of functions {0„}in D, which we say they go to 

zero in D, and write it as 0„— > 0. For any such sequence, we require that it satis- 
fies two conditions: 

i- There is a bounded interval J such that the support of 0 „ c J for all n, and 

ii- 0 and all its derivatives go to zero uniformly in J . 

We make two useful comments here concerning these conditions. First, in general, 
different sequences going to zero in D have different bounded interval J associated 
with them. In other words, we require only that for a given sequence going to zero 
in D such a bounded interval to exist. Second, condition ii simply means that if m is 
the order of the derivative, m - 0, 1,2, ...., and given s > 0, we can find N(m) 
such that for any x in the interval J, we have |0^ m) (x)| < s for all n > N(m). Note 
that the word “uniformly” means that the number N(m) is independent of x but it 
can depend on the order of differentiation m. Here we have used the notation 0 ( ' M) (x) 
for the mth derivative of 0„(x). 

Example 1- We can show easily that the following sequence goes to zero in D: 
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l 

0„(x) = - 0 (x; a) (See eq. (4) for the definition of 0 (x; a)) (6) 

n 

Example 2- The following sequence of functions does not go to zero in D: 

1 (x \ 

0 n( x ) = ~ 0 - ; (See eq. (4) for the definition of <p (x; a)) (7) 

n \n ) 

The reason is that for a given function <p n , the support is the interval [-na, na ] 

and, therefore, the support cannot remain bounded for all n. Condition ii is, thus, 

violated. 

Now we can give the definition of continuity of a functional: 

D 

A functional F[cf)] is continuous if for any 0 n — > 0, we have F[(p n ] — > 0 . 

Note that { F[(p n ] } is simply a sequence of real or complex numbers. By using the 

properties of (Lebesgue) integrals, we can show the following very important result: 

Any OF /(x) generates a continuous linear 
functionlthmxxgh the relation : F[(p] - J)^ f(x) (f) (x) d x 
(0 is in the test function space D ) 




F[(j)] =/f(x) (j)(x) dx 

Figure 2- The functional rule used to define an OF /(x) by a table of its func- 
tional values over the space D 
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Figure 2 shows schematically how the table of functional values is constructed for 
an OF. So we can use the above functional (or rule) to construct the sequence or the 
table of numbers {F[f], (fie D\ to identify (or describe , distinguish ) the OF fix ) . 
Are there any more continuous linear functionals described by other rules? The 
answer is affirmative and, in fact, there are many of them. Let us consider the follow- 
ing example: 

Example 3- Define the functional 8 [<fi] = 0 (0), 0 eD. We can easily show that this 
functional is linear. We now show that it is also continuous. Take any sequence that 

goes to zero in D, say {(p n ). Since 0„— > 0, then in particular <p n — >0 uniformly. The 
uniformity implies that it does not matter where x is. Therefore, 0„(O) = 6 [f n ] — >0. 
Can you see how we use condition ii of the definition of continuity of a functional 
here? Because of the importance of this functional in GF theory, we will give it the 
name of the Dirac functional. This is obviously related to the Dirac delta function. 
We will establish the connection below (see eq. (9)). 

Using elementary Lebesgue integration theory, we can show that there is no OF 
d(x) that can generate the continuous linear functional 6 [<fi] by using the rule given 
in eq. (3). Therefore not all the continuous linear functionals are generated by OFs 
based on eq. (3)! 

Now if we decide to define (or describe, identify) any OF conceptually by the table 
of its functional values using the rule of eq. (3), we may just as well say that the 
table generated by any continuous linear functional also defines a function. We 
now have new objects (tables) that cannot be identified with OFs, e.g., the func- 
tional in Example 3. We call the entire collection of objects (tables of functional 
values) generalized functions (GFs). We have achieved our goal of generalizing the 
definition of OFs. We, thus, state the following: 

Generalized functions are defined by the tables 
of functional values of continuous linear functionals. 

Since we are now going to think of functions in the new way by their functional 
values, we will use a slight abuse of language and say: 


Generalized functions are continuous linear functionals . 
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The particular space of GFs that we have constructed using the test function space 
D is known as the (Schwartz) GF space D' . We know that by construction, all OFs 
are also GFs. We refer to OFs as regular generalized functions. The new objects 
(functions) that cannot be identified as OFs are called singular generalized func- 
tions (Think of the GF in Example 3 now.). Figure 3 shows how one should think 
of the space of GFs. In applications, several other test function spaces are used but 
the idea behind the construction of the associated GF spaces is exactly the same as 
we have presented above. At this stage, it is important to know that there are many 
singular GFs that are absolutely essential in applications 7 . As a matter of fact, the 
introduction of the singular GFs into mathematics and the development of a rigor- 
ous operational algebra to work with them have provided a very powerful analytical 
tool to physicists and engineers. 



Singular GFs 


be D 


G[<(>] =j g(x) p(x) dx 
g(x) symbolic function 


Real or complex 
numbers 


Figure 3- The space D' of GFs consisting of regular GFs (OFs) and singular 
GF that are new 


Readers should be aware that the definition of GFs given above, like the construc- 
tion of real numbers by Cauchy convergent sequences in analysis, is of enormous 
conceptual value. We will have to use this definition to obtain the properties of GFs 
and to extend the useful operations on OFs to all GFs. But once we have achieved 
these goals, we want to get away from dealing with functionals as functions and 
work with GFs with the same notations and ease as OFs. Like the man in the street 
who never needs to know about convergent sequences to work with numbers, we 
want to use GFs without being burdened by the conceptual framework of the con- 
struction of GFs. We can actually do this but first let us look at an example of a GF. 

Example 4- Here is a continuous linear functional defining a GF: 
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AX) 

G[0] - I e 2 x sin 3 x 0 (x) dx - 5 0 (0) 0 e D (8) 

J—oo 

You should have little problem proving this fact. Note that, in practice, we always 
give the rule (the definition of the functional) only and we do not attempt to con- 
struct the table of the functional values because it is technically impossible. But we 
can use the rule to make part of the table if we want to ! 

We will now dispose of the functional notation by the introduction of symbolic 
functions. This is a brilliant idea. It will allow us to work with singular generalized 
functions just like the OFs. We start with the Dirac functional of Example 3 which 
is a singular GF. We introduce a symbolic function <5(x) which is supposed to hold 
the information (or the memory) about the Dirac functional when we are doing 
algebraic manipulations until it appears as an integrand of an integral in a product 
with a test function .We then use the following relation agreed upon by convention : 



0 (x) d(x) dx = d[0] = 0 (0) 


0eD 


(9) 


It is very important to recognize that ‘the integral’ on the left of this equation is not 
a Riemann or a Lebesgue integral because d(x) is not an OF. You should think of 
the whole integral symbol like a Chinese character standing for d[0] = 0 (0) which 
is now completely meaningful. The symbolic function d(x), known universally as 
the Dirac delta function, can be used in algebraic manipulations like an OF until it 
appears under an integral sign like that on the left of eq. (9). You should, however, 
learn some more about the properties of singular GFs before you can comfortably 
or boldly manipulate them algebraically. We will give you what you need to know 
later but first let us look again at Example 4 using symbolic function notation. 

Example 5- The GF in Example 4 can be written symbolically as 
g(x) - e 2x sin 3 x - 5 d(x). This GF is the sum of a regular and a singular GF. 


2.2 - Multidimensional GFs 

The above definition of GFs can easily extended to multidimensions. The test 
function space D is now n-dimensional and consists of all infinitely 
differentiable functions with bounded support. An example of such a function 
can be constructed from the function in eq. (4). Let x = (xi , X 2 , x n ), and 

I x\ = (x? + x 2 + ... + x^) 1,/- . Then the we can show that for a > 0, the 
following function is in space D: 
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( 10 ) 


From this function and any continuous function, we can generate another test func- 
tion by convolution over a finite region of space. Therefore, the space of test func- 
tions D has an uncountably many functions in it. We then define the multidimen- 
sional (Schwartz) GF space D' as the space of continuous linear functionals on 
D. This definition is appropriate for discovering the operational properties of GFs 
that give us a powerful tool for applications. In practice, however, we will work 
with symbolic functions when we deal with singular generalized functions. There- 
fore, we will be always working with generalized functions as if we are working 
with OFs. Intuitively, this is of great help to people interested in applications. The 
situation is similar to the use of identical notations for real and complex numbers as 
much as possible when we are solving a problem in the complex plane. 

From our point of view, the multidimensional generalized functions are much more 
important than GFs of one variable in studying partial differential equations, particu- 
larly in wave propagation problems. It is important to recognize that to study the 
multidimensional GFs, one needs to have a working knowledge of differential geom- 
etry of curves and surfaces 30 . 

As in the case of one variable, the most important singular GFs in multidimensions 
are the Dirac delta function and its derivatives. Of interest to us is the delta function 
with the support on a surface which can be in motion. This function has no one 
dimensional analogue. We start with a simple example and then discus the Dirac 
delta function with the support on a surface. 

Example 6- The Dirac delta function b(jr) with support at the origin has the follow- 
ing property 


where V is an arbitrary volume which includes the origin. 

Example 7- Referring to Figure 4, let S : f(x) = 0 be a smooth surface defined in 
such a way that the unit outward normal n = V /. This can always be done because 
if n + V /, then define this surface by the new implicit function /( x) / 1 V / | = 0 
which has the desired property. 



( 11 ) 
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Figure 4- The description of the surface S defined implicitly by f(x) - 0 with 
unit outward normal n = v/. 


We want to interpret the important GF 6(f) often appearing in applications. We 
define a Gaussian curvilinear coordinates (u l , it 2 ) over the surface S and define the 

third coordinate u 3 - f along the local normal to S to all points in the space in the 
vicinity of S. Note that we have 


d x- 


V 8(2 )(M) 


d u 


u 


= (w 1 , 


u 2 , u 3 ) 


( 12 ) 


where g( 2 )(u) is the determinant of the first fundamental form of the surface 
u 3 = /( x) - constant 29-31 . We thus have the following very useful result 


fm ^ w d* - fa) 0 wo facin' d * d j = 

J< <p u 2 , 0)) yj g[ 2 ){u l , u 2 , 0) du l du 2 - J" <p(x)dS 


(13) 


You can see the importance of differential geometry here. Note that we have used 
the definition of the Dirac delta function in one variable when we integrate with 
respect to variable u 3 after the second equality sign. For a more geometric deriva- 
tion of this result see the papers by Farassat 29 ’ 31 . 
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It is important to remember that one cannot define a GF at a single point but only 
over a finite interval. This does not cause any difficulties in applications. We can 
talk about equality of two GFs over an interval. We say that two GFs F[f] and G[0] 
are equal over an interval I if for all 0 e D with support fully within I, we have 
F[f] = G[(p]. Remember, a functional is a rule that gives us a number for any func- 
tion in D. Therefore, we are talking about equality of two numbers obtained from 
two rules. From this definition we can say that 8 (x) = 0 (the zero function and not 
the real number zero) over the interval (-oo, -s] (J [s, oo) for any s > 0. 

To learn more about the multidimensional GFs we recommend the books by 
Gelfand and Shilov 7 , Stakgold 8 , Vladimirov 9 , Kanwal 11 , Jones 12 , Fenyo and Frey 19 , 
and the two NASA papers by Farassat 29,31 . For more advanced treatment as well 
some very important recent advances, see the books by Hormander 5 and Taylor 6 . 


2.3- Operations on GFs 

Now that we have extended the definition of OFs to GFs, we must extend the defini- 
tions of operations on OFs to GFs. The general approach is as follows. Write the 
definition of the operation in the language of linear functionals and then use it for 
all GFs if it makes any sense (see the definition of generalized Fourier Transform 
below where we have to change the test function space). In some cases the exten- 
sion of the definition of an operation gives us the opportunity to solve new prob- 
lems or solve old problems in much simpler ways. We single out generalized differ- 
entiation here as an example of an operation of considerable importance in applied 
mathematics. 

In the following we assume that all test functions are in space D. 


2.3.1- Simple Operations 


1- Addition of two GFs- If f(x) and g(x) are two OFs described by functionals 
F[f \ and G[0], 0 e D then using the linear property of integrals, we can write 




Jffdx+ J'i 


(f + g)4>dx - I / 0 dx + g 0 dx = F[0] + G[0] 


(14) 


Basically, we know the left side of the first equality is the functional representation 
of the sum of the OFs f(x) and g(x) and is given by the right side of the last equality 
sign. The above result now makes sense for all GFs, i.e., we define 


(F + G)[0] = F[0] + G[0] 


( 15 ) 


Note that (F + G) is a single functional which stands for the sum of the two GFs. 
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Written in symbolic function notation, there is no difference between the notation 
for OFs and GFs for the sum of two functions. For example, for the sum of an OF 
fix) and the Dirac delta function we can simply write fix) + Six). 

2- Multiplication of two GFs- This operation cannot be extended to all of 
Schwartz GFs. In fact it is the aim of the new nonlinear GF theory to be able to 
define multiplication for all GFs 22 ’ 23 . Any two OFs can be multiplied. We can 
define multiplication of any GF with any C°° function. Again let fix) be any OF 
and let a(x) e C°°, then 


J' a(x) f(x ) f (x) 



(a(x) fix)) dx = F[a f] 


(16) 


We know that the left side of the first equality sign is the functional representation 
of the OF aix) fix). Since the function a feD, the expression on the right of the 
third equality sign is a functional on the test function space D. The rule of multiplica 
tion of a generalized function F[f] and a C°° function aix) is given by the rule 


aF[f ] = F[a f] 


(17) 


Note that in this equation, the symbol aF does not mean the product of a and F but 
stands for a single functional, i.e., a rule described by the right side of eq. (17). 
Again in symbolic notation we simply write aix) fix) but now we have a rule to 
interpret such products when fix) is a singular GF. We give an example here. 

Example 8- What is aix) d(x) where a e C°°? Let f e D, then 


f 


a <5[0] = I aix) fix) d(x) dx - 5 [a f] = g(0) f (0) = I «(0) fix) d(x) dx 


!‘ 


(18) 


We can write this result in symbolic function notation as 


aix) Six) = ai 0) d(x) 


(19) 


This useful result has been known for many years to physicists. 

3- Translation of a GF- Let fix) be an OF, and let us define E h fix) - fix - h), 
i.e., translating the function to the right for h > 0. In functional notation, we have 




f 


E h fix) fix) dx = fix- h) fix) dx = fiy) fiy + h)dy = F[E_ h f] 


5 > 


( 20 ) 


Again the left side of the first equality sign is the functional generated by E h fix) 
for which we use the symbol E h F[f\. Since for feD, we have E_ h feD. Therefore, 
we use the following result for the definition of the translation of any GF: 

E h F[f] - F[E_ h f] 


( 21 ) 
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Once again we have obtained a nontrivial result. In symbolic function notation we 
use the same notation for GFs as for OFs but occasionally we have to resort to eq. 
(21) to interpret the exact meaning of the translated GF such as in the following 
example. 

Example 9- What is 5(yc - a) where a is a real constant? The functional representa- 
tion of this GF is obviously E a 5[<p]. Using the above rule, we have 


f 


6(x - a ) 0(x) dx - E a d[0] = S[E_ a </>] = 5[<fi(x + a)] = (p{a) 


( 22 ) 


Again this result has been known by physicists for a long time. 

We can extend the definition of other simple operations on OFs to GFs such as the 
expansion and contraction of the scale of the independent variable. For example, 
one can rigorously show that for a ± 0 we have the following useful result 


1 

6(a x) = — d(x) 

\a\ 


(23) 


2.3.2- More Advanced Operations 


We will extend the definition of two operations on OFs to GFs here. These are 
Fourier transformation and differentiation operations. 

1- Fourier Transform of GFs- We work in one dimension here. Let f{x ) be a func- 
tion defined on R which has Fourier transform defined by the relation 

/O') = r f(x)e 2 ’ ixy dx (24) 

kJ—OO 

If the functional F[<p] is identified with f(x ) , then its Fourier transform must be 

A 

identified with f(y ) as follows 


^[0] = f(y) 0 00 dy = f{y ) 0 (y) dy = F[0] 


(25) 
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where 0 is the Fourier transform of 0. We have used Parseval’s theorem after the 

/V /V 

second equality sign. Now the functional F[f] is only defined if 0eD whenever 

A 

0 eD. We can show that, in general, this is not so and 0 can lie outside the space D. 
We need to define a test function space S on R consisting of C°°functions which go 
to zero at infinity faster than \x~ n \ for any n. This can be done and in fact we can 
show that D c S. The space of tempered GFs S' is the space of all continuous lin- 
ear functionals on S. An important relation to remember is D c S c S' c D' . All the 
GFs in the space S' have generalized Fourier transform given by the relation 

y\ /v 

F[0] = F[f] where now feS. We will say no more on this interesting subject here. 
See the references at the end of this paper. We strongly recommend the book by 
Strichartz 10 for this subject. 

2- Differentiation of GFs- This is the most important operation on GFs because of 
its many applications in mathematics. Many of the problems of wave propagation 
need the use of this concept. Let us start with a C 1 function f(x) which we identify 
with the functional F[f ] . The functional representation of f'ix) can be manipulated 
as follows 


F'[(f>] 



\x) 0 (x) dx 



-F[<f>'] 


(26) 


Note that we have integrated by parts after the second equality sign and we have 
used the fact that the support of 0 is compact to drop the terms involving the limits 
of the integral. Since we know that the function 0' is in D, F[<p'] is a proper func- 
tional on D and we can use the following relation for the definition of the general- 
ized derivative (GD) of a GF 


F'W = ~F[0 1 

From this we can define the nth generalized derivative of a GF 

F (n) [0] = (-1 ) n F[f {n) ] 


(27) 

(28) 


This significant result states that GFs have GDs of all orders. When we work with 
symbolic functions we often utilize a bar over a derivative operator to designate 
generalized differentiation if the operation can be confused with ordinary differen- 
tiation. See eqs. (30) and (34) below. For more on this notation see Note 2 at the 
end of this paper. It is understood that the differentiation of singular GFs can only 
be GD and in that case we will not use this notation. 


Example 10- GD of the Heaviside function h(x) - Oforx < 0 and 
h(x ) = 1 forx > 0 is found as follows: 
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AX3 

H'[(f)] = -H[(p f ] = - J (f)' (x) d x - <p (0) 


(29) 


Note that we have used the fact that the test function has a bounded support and, 
thus, we have <p(oo) = 0. Therefore, we have 


h\x ) = <5(x) 


(30) 


That is, t/n? GD o/t/ze Heaviside function is the Dirac delta function . Note that 
h\x) - 0 but h\x ) ^ 0. We will return to this point later. This result shows that the 
GD of an OF can be a singular GF! 

Example 11- The generalized derivative of the Dirac delta function is 


6'[f] = JS'(x) (f> (x) dx 
Similarly, we have 
6 (n) [(f)] = Js 0i) (x) ([) (x) 


= -5W] = -0'(O) 


= (-l) n 0 (n) (O) 


(31) 


(32) 


Note that in eqs. (31) and (32), the integrals are meaningless; they simply stand for 
the functional rules to the left of the first equality sign in the respective equations. 

Example 12- Let the OF f(x) be differentiable with a single jump of 
A / = /(xo+) - f(x o-) at the point xo as shown in Figure 5. 



Xo 


Figure 5- A differentiable function with a single jump discontinuity at xo 
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We find the GD of this function as follows. Let F[f ] be the functional representa- 
tion of this function. Then its GD can be manipulated as follows 

F'[f\ = 


f 


-F[<f>'] = - fix) (f>'(x) dx = - 


r 

J—oo 


r^oo 

f(x ) fix) dx - I f{x) fix) dx = 

Jx 0+ 


rx o- 

J r 


(x) fix) dx + fix) fix) dx + (/(x 0+ ) - fixo-)) fx 0 ) = 
Jx 0 + 


(33) 



(x) fix) d X + A / fix o) 



if'ix) + A / d(x - xo)) f (x) dx 


Interpreted in symbolic function notation, we have obtained the following important 
result 


fix) = fix) + A / Six - x 0 ) 


(34) 


Again we have shown that the GD of an OF can be a singular GF. Furthermore, 
this result shows the need for a notation for generalized differentiation which we 
have utilized here, i.e., the overbar on the left side. 

Example 13- Referring to Figure 4, let the function kix) be a differentiable func- 
tion with a jump discontinuity across the surface f - 0. We define the jump of the 
function kix) across the surface as follows: 

A k(u l , u 2 ) - k(u l , u 2 , 0+) - k(u l , u 2 , 0_) = k(x |/=o+) - k(x |/=o-) (35) 

where u - (V, u 2 , n 3 ) is the curvilinear coordinate system we used in Example 7. 

Note that the function kix) is discontinuous in variable u 3 - f only. This is pre- 
cisely why we selected such a coordinate system in the vicinity of the surface of 
discontinuity. Now we are going to use eq. (34) to derive the following important 
result: 


dkixiu)) 

dxj 

df dk 


dx, du j 


du 3 , ox dk df dk 

+ A k 6 ( u 3 ) = — f A k — — 6 if) - — h A k n\ 8 if) 


(36) 


dXr 


dx-. 


dx-. 


dx. 


We have used the summation convention on repeated index in the above equation. 
We can write the above result in vector notation as follows: 


V k (x) = V kix) + A knd if) 


(Note n = V/ here.) 


(37) 
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We give here two more important relations which can be remembered easily. Let 
k (x) be a vector field with a jump discontinuity across the surface f - 0 with the 
jump defined as in eq. (35). Then the following results hold: 


v • k (x) = v • k (x) + A k ■ n 5 (/) 


(Note w = V/ here.) 


(38) 


v x k (x) = v x k (x) + n x A k 8 (/) 


(Note w = v/here.) 


(39) 


See Note 3 for a comment on Eqs. (37) to (39). With what we have obtained so far 
we can solve a number of interesting problems. Before we give some examples, let 
us say a little more about why one should use the generalized derivative. 


2.4- Why We Use Generalized Derivatives 

Part 2 of this tutorial is on applications of generalized derivatives and we will have 
much to say about GDs. Here we mention some useful results that readers should 
know. The most important fact to remember about GDs is that in GF theory gener- 
alized differentiation is a continuous operation. In particular, generalized differen 
tiation and integration are inverse operations. This means that the integration of GD 
of a GF recovers the function faithfully with all of its jumps. It is well known that 
this is not so in the case of the ordinary derivative of an OF. 

Fet /(x) be a discontinuous function as in Example 12. The symbolic meaning of 
J* /' (y) dy can be found as follows. Form the functional 

HILT 00 dy) 000 dx =JT r 00 { X°°0U) dx} dy. Now define 

if/(y) = jT° </>(x) dx which is in the test function space D and clearly we have 

if/' (x) = -(f) (x). From Eq. (33) we have 

tj’ O’) '/'O’J dy = - Jy /(A-) '// (x) dx = C jtx) <l>(x) dx. Symbolically, this means 

that f„ /' (y) dy = /(x), even if the function /(x) is discontinuous as assumed. Note 

that we have essentially followed the procedure that Schwartz used to define the 

indefinite integral of a generalized function 4, 7 . 

Example 14- Fet h(x) be the Heaviside function defined in Example 10. Then from 
the above result we conclude that 



h ' (y) dy 



6(y) dy - h (x) 


but 



h ' (y) dy - 0 =£ h (x) 


(40) 
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Therefore, any problem involving unknowns with discontinuities should be set 
up in GF space and all derivatives must be treated as GDs from the start. Such 
discontinuities are either artificially introduced, e.g., as in FW-H equation, or are 
genuine, e.g., a shock wave in fluid flow field. The statement following eq. (40) has 
profound implications for applications some of which will mentioned below: 

1- All local conservation laws obtained from control volume analysis are valid 
if all the derivatives are considered GDs. Examples of such conservation laws 
are the mass continuity and the Navier-Stokes equations. 

2- The jump conditions across a natural discontinuity in a flow field are inher- 
ent in the conservation laws. These jump conditions are obtained easily by apply- 
ing the definition of GD to the conservation laws. See Example 15 below. 

3- Any problem involving derivatives in which singular GFs appear must be 
set up in GF space. One must then carefully distinguish between the GDs and the 
ordinary derivatives of any function appearing in the algebra. An example of such a 
problem is finding the Green’s function of an ODE or a PDE. 

4- The generalized differentiation operator commutes with all common limit 
operations. This is the most pleasant and useful property of generalized differentia- 
tion. The following operations are allowed provided that the limits of the integra- 
tion are independent of the variable of differentiation: 


d 


dxj 



q (x, y ) dy = 



dq(x,y) 

dxi 



lim q n {x, y) = lim 

n^oo n->oo 


dg n ix, y) 
dxj 


(41) 


See the NASA paper by Farassat 31 for more examples. 

Example 15- Let us find one of the jump conditions across an unsteady shock wave 
which is described by the surface /( x, t) - 0. We define the jumps in fluid parame- 
ters as in eq. (35). The mass continuity equation is 

dp 

— + v • (p u) = 0 (42) 

dt 

Here the fluid density is p and the fluid velocity is u. First we know this conserva- 
tion law was derived from a control volume analysis. Therefore, it is valid if we 
treat all ordinary derivatives as GDs, i.e., we have 

dp _ 

— + V • (p u) = 0 
dt 

We now use the definition of GD in this equation as follows 


( 43 ) 
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dp 

dt 


+ v • (p u) = 


dp df 

— + A P — £(/) + v • (p w) + A(pu)-nd(f) = A[p (u n 
dt dt 


Vn)]d(f)= 0 


where we have used the notations u n - u n and v„ = -df /dt is the local normal 
velocity of the shock. We have used eq. (42) in the above result. Therefore, one 
jump condition across an unsteady shock wave is 

A[p (u n - v„)] = 0 or A(p w„) - v„ A p = 0 (45) 


Note that we derived this condition without using the pill box analysis employed in 
the classical derivation of the result. 


5- We can use the Green’s function technique to find the discontinuous 
solutions of an ODE or a PDE provided that we formulate the problem in GF 
space and use generalized differentiation everywhere. This is a very powerful 
result and its impact on applications is enormous. We will not give any example 
here but you can see some good examples in the next section. 

We conclude this section by stating the fundamental theorem that characterizes 
GFs in D' : 


6- GFs in D' are GDs of finite orders of continuous functions 7 . This is a very 
profound result obtained by Laurent Schwartz. As a simple example we note that 
the Dirac delta function is the generalized second derivative of the continuous func- 
tion that is zero on the negative axis and is equal to x on the positive axis. 

We will elaborate on all these results in part 2 of the present tutorial. 
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■ 2.5- A Note on Test Function Spaces 

We have developed above the GF theory using the test function space D. Although 
the GF spaces D' and S' are very useful in applications, we often have to use test 
functions other than the spaces D and S which are not C°°functions. But the con- 
struction of GFs in the space D'can be used as a model to get other kinds of 
GFs. For example if we take a space of test functions consisting of C n functions, 
then we can define GFs on this space as continuous linear functionals which only 
have GDs up to the order n. We must, of course, modify our definition of the conti- 
nuity of a linear functional appropriately. To define GFs suitable to treat ODEs and 
PDEs, say a differential operator L (a differential equation plus linear homogeneous 
boundary conditions (BCs)), the test functions must have continuous derivatives of 
some order (related to the order of highest derivative of L) and should satisfy the 
adjoint BCs. All things we have learned about GFs here, particularly the operations 
on GFs, apply as they are. For this reason we do not specifically identify the test 
function space unless such identification will shed some light on the problem at 

hand. From this discussion, it is obvious that we can write J)^6(x) dx - 1 even if 

the test function f(x) - 1 in the integrand is not in the test function space D. See 
also Example 14. 

3- Some Applications of Generalized Differentiation 

We will give two important applications of generalized differentiation which are 
very useful in solving wave propagation problems. These are: 

1- The Green’s function of a second order ODE, and 

2- The imbedding of a problem in another problem whose Green’s function is 
known. 

You should note that in these applications we have to utilize almost everything we 
presented on GFs in the previous section. 

■ 3.1- The Green’s Function of a Second Order ODE 

In working with wave propagation problems in the time or frequency domains, we 
often have to obtain the Green’s function of a second order ODE in intermediate 
steps. Since we are now able to explain some of the subtleties of the derivation of 
this Green’s function from the point of GF theory, we will present a discussion of it 
here. We hope that you will agree with us that GF theory is the right tool to use 
here. 
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Consider the following second order ODE with two linear homogeneous boundary 
conditions (BCs): 

/ u = a(x) u"(x ) + b(x) u'(x ) + c(x) u(x) = k(x) x e [0, 1] (46) 

BCi [w] = 0, BC 2 [u] = 0 (47) 


The Green’s function g(x, y) for this problem is formulated classically as follows 8 
(see volume 1 of this reference): 


l x g(x, y ) = a(x) 


d 2 g(x, y) 

dx 2 


+ b(x) 


dg(x, y) 
dx 


+ c(x)g(x,y ) = d(x-y), 


(48) 


x, ye [0, 1] 

BCi[g(x, y)] = 0, BC 2 [g(x, y)] = 0 (in variable x only) (49) 


Now we see that we have the Dirac delta function, which is a singular GF, on the 
right of eq. (48). So the partial derivatives in this equation are most likely GDs (we 
know the answer but we are presenting what most books for engineers and physi- 
cists give). But how did we suddenly get into the GF space? Well, we did not really 
set up our problem correctly. First from the theory of ODEs we must know some- 
thing about the nature of the solution of the BC problem described by eqs.(46) and 
(47) . Fet us assume that we know that ueC 2 . This can be proved a posteriori when 
the solution is at hand. In that case, in the domain of the definition of this function, 
since a C 2 function cannot have any jumps or sharp corners, we have 

/ u = a(x) Ti"(x) + b(x) 1 f (x) + c(x) ii(x) = ^ 

a(x) Li'(x) + b(x) u'(x) + c(x) u(x) = l u 

The Green’s function is used to find the unknown function from the following 
relation: 


u(x) = f k(y) g(x, y) dy 
Jo 


(51) 


Fet us discover further properties of the Green’s function from this relation by apply- 
ing the linear differential operator (the ODE plus the BCs) to both sides of the equal- 
ity. First we get 


/ u(x) = l Ii(x) = l x I k(y) g(x, 
Jo 


y)dy = f k(y) l x g(x, y) dy = k(x) 
Jo 


(52) 


Note that we freely took the generalized derivatives in l x inside the integral sign 
without concern because this step is allowed in GF theory. From this we immedi- 
ately can see that 
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lx g(x, y ) = a(x) 


d 2 g(x, y ) 
dx 2 


+ b{x ) 


dg(x,y) 

dx 


+ c(x) g(x, y) = 6 (x - y) 


(53) 


Therefore, we established rigorously that all the derivatives in eq. (48) are GDs. 
We will get back to the implications of this fact soon. Next we apply the two BCs 
to both sides of eq. (51) remembering that these are linear and homogeneous: 

BC V [w(x)] = BC X f k(y) g(x, y) dy = 

Jo 

where the symbol BC A stands for either of the two BCs but applying to variable x 
only. From the above result, because the function k(y) is arbitrary, we conclude that 
the Green’s function in variable x satisfies both BCs: 


f 


k(y ) BC V [ g(x, y)] dy = 0 


(54) 


BCi ; V [g(x,y)] and BC 2>x [g(x, y)] 


(55) 


Now let us see what eq. (53) implies about the Green’s function. Obviously this 
function changes its character at x = y. So let us assume that it is described as fol- 
lows by two C 2 functions with possible jumps in their values and their derivatives at 

x = y- 


g(x, y) 


(gi(x,y) 

\ gi(x, y) 


x<y 

x>y 


(56) 


First of all since the Dirac delta function 8 (x - y) is identically zero in the two open 
intervals on both sides of x = y, we must have: 


l x gi(x,y) = l x giix,y) = 0 


(57) 


This means that the functions gi(x, y) and g 2 (x, y) in variable x are solutions of the 
homogeneous ODE l x g(x, y) - 0. Let us take the GD of the first and second order 
of the Green’s function described by eq. (56) 


dg(x, y) 
dx 

d 2 g(x, y) 
dx 2 


dg(x, y) 
dx 


+ [g 2 (y+, y) - gi(y~, y)] S(x-y) 


d 2 g(x, y) 
dx 2 


dg 2 (y+,y ) 
dx 


dgi(y~,y) - 

dx 


8(x-y) + Lg 2 ()4, y) - gliy-, y)] s' (x - y) 


(58) 

(59) 


We must substitute these two relations in eq. (53) and equate both sides term by 
term. Since there is no derivative of the Dirac delta function on the right of eq. (53), 
we immediately conclude that 
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#2(;y+,;y)-giCy-,;y) = o 


(60) 


This means that the Green’s function is continuous at x = y. Furthermore, we get 


rd^O^) 

<9gi(y_,y)i 

1 

L dx 

dx 

a(y) 


( 61 ) 


(assuming that a(y) ± 0forye[0, 1]) 

This means that the derivative of the Green’s function with respect to x has a jump 
described by the above equation at x = y. We can show that eqs. (55), (57), (60) and 
(61) uniquely define the Green’s function. Note that only by using GF theory have 
we been able to get the above results rigorously and satisfactorily. 

Example 16- Let us take l u(x) = u " for x e [0, 1] with the following BCs: 

u (0) — 2 u' (O') — 0 and w(l) + w'(l) = 0 (62) 


The above procedure will give the following Green’s function for this problem: 

(y ~2)(x + 2) 

g(x, y ) = 


x < y 


Cy + 2) (x - 2) 


(63) 


X > y 


This Green’s function satisfies the conditions derived above. Note that the 
differential operator in this example is selfadjoint and, therefore, the symmetry of 
the Green’s function in variables x and y is expected 8 . 


3.2- The Imbedding of a Problem in Another Problem 

We became interested in this problem through our work on the Ffowcs Williams- 
Hawkings (FW-H) equation 2 . The situation is described in Figure 6 below. 
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this data surface 

Figure 6- The problem of noise radiation from a moving body (a helicopter 
rotor blade here) as formulated by Ffowcs Williams and Hawkings 


We are given an inhomogeneous wave equation in the exterior of a moving surface 
described implicitly by the function /( x, t) = 0. This surface will be called the data 
surface here. Ffowcs Williams and Hawkings searched for a radiation formula in 
the exterior of the data surface. Essentially, these authors used an imbedding 
method by making their wave equation valid in the entire unbounded space 2 . This 
imbedding allowed the use of the free space Green’s function for getting the radia- 
tion formula. They assumed that the fluid inside the data surface is always at rest 
with the condition of the undisturbed medium. Since an artificial discontinuity has 
been introduced in the medium, one must work with GFs. This imbedding tech- 
nique provides us with an elegant method of deriving a radiation formula. This 
technique can be applied to a whole range of other problems involving ODEs and 
PDEs, particularly those related to wave propagation. We will discuss some of 
these problems here. First, however, we present the imbedding technique. 
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3.2.1- The Imbedding Technique 

We will consider the problem of radiation from a moving and deformable surface 
/( x, t ) = 0, assuming that V / = n, where n is the unit outward normal. We have 
shown the moving surface at two times in space in Figure 7. 



Figure 7- A moving deformable surface described by f(x, t) = 0 at two differ- 
ent times ?iand t 2 in space 

We want to solve the following exterior acoustic radiation problem: 

1 d2(f V72 

V = 

c 1 2 dt 2 (64) 

n 2 (p - Q (x, t ) (jt in exterior of / = 0, c = constant) 

We are assuming that we have all the data about the function on the moving 
surface. These data may be overdetermined but right now we are not concerned 
about this fact. We want to extend the domain of the problem to the interior of the 
surface so that we can use the simple free space Green’s function of the wave 
equation. In other words, we want to imbed our problem in another problem with a 
larger domain in space (this can be done also for the dimension of time). This 
imbedding will require the powerful machinery of GFs. Here is how we do the 
imbedding. We start by extending the definition of the two functions in eq. (64) as 
follows: 
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P (x, t) = 


(65) 


<pix,t) f > 0 ~ ( Q(x,t) f > 0 

, and Q (x, t ) = { 

0 /< 0 (o /<0 

The first thing we notice is that with this extension, the new functions satisfy the 
following wave equation: 

□ 2 p = Q (x, t) (66) 

Note that the unknown function p is discontinuous. To find this function, we must 
set the problem in GF space. This simply means that we find what the wave opera- 
tor looks like when all the ordinary derivatives are written as GDs as follows: 


dp dp) df dp 

— = — + — 5(f) = (pv n 6(f) 

dt dt dt dt 


(67) 


where v n = - df /dt is local normal velocity of the surface which is unambiguously 
defined even for a deformable surface. Also we use the notation ip - ip (f - 0+) 
when this function multiplies the Dirac delta function 6 (f) as in the last term on the 
right of eq. (67). Next we find the following relations: 


d 2 p d 2 p dip d 

— v = — y - — v » 5(f) - — [ip v n 6(f)] 
dt 2 dt 2 dt dt 


( 68 ) 


v 2 p = v 2 p + n-Vtp 6(f) + v • [ip n 6(f)] (69) 

where we have again used the notation that any function of p multiplying the Dirac 
delta function is evaluated at / = 0+ . Let us define the local normal Mach number 
of the surface as M n = v n /c, and define additionally the following two symbols: 
p> t - dip I dt and p> n = « • V ip. Then the above two relations will give us: 


n 2 p = u 2 p- - ip t M n 6(f) - - [ip M n 6(f)] - ip n 6(f) -v -[ipn 6(f)] 
c c dt 

Or, after using eq. (66), we obtain 


(70) 


( 71 ) 




FF&MKM- GFs Pt 1 (11-09) V5.nb 


34 


Note that we do need to use a bar over the wave operator on the left side since from 
the right side it is clear that all derivatives involved are GDs. Now we can use the 
Green’s function of the wave equation in the unbounded space, the so-called free- 
space Green’s function, to find the unknown function <p(x, t) everywhere in space. 
The result is the Kirchhoff formula for moving surfaces 31,32 . The classical deriva- 
tion of this formula by Morgans 33 is very complicated and in fact until Farassat and 
Myers rederived it by the modern method presented here, there were some doubts 
expressed about the correctness of Morgans’ result 2,34 . 

We will not derive the Kirchhoff formula for moving surfaces here since the two 
references by Farassat and Myers 32 and Farassat 31 are quite comprehensive alleviat- 
ing the need for further elaboration. However, we next derive some classical results 
for stationary surfaces by the imbedding method which are surprisingly simpler 
than the classical methods of their derivation. 

Before we end this discussion we want to mention an important question that comes 
to mind. Answering this question may lead to further applications of the Green’s 
function of the unbounded space. To imbed our exterior radiation problem in R 3 we 
have imposed that inside the surface both f and Q be zero. Can the functions f 

and <2 be assumed to take other values and are there advantages in doing so? 
We will find the answer to this question using an ODE to simplify the analysis. 

■ 3.2.2- Derivation of Some Classical Results Using the Imbedding 
Technique 

1- The Green’s Theorem for the Laplace Equation in the Exterior Domain of a 
Surface- The problem is described in Figure 8. 
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Obs. 


Figure 8- Problem description for deriving the Green’s theorem for Laplace 
equation, / > 0 outside the surface and n - v / on the surface 

We imbed this problem in R 3 by extending the definition of <p to the interior of the 
surface as in eq. (65). We find generalized Laplace equation next: 

v 2 <£ = 

- ( 12 ) 

V 2 (p + n- V p 8(f) + V -[(pit 6(f)] = (p n 8(f) + v-[(pn 6(f)] (xeR 3 ) 

The definition of all symbols are identical to those in Section 3.2.1. We have used 
the fact that V 2 f = 0 in R 3 . The Green’s function of this problem in the unbounded 
space is - 1 / (4 n r), where r - \ x - y \ . Here x and y are the observer and the source 
(field) variables, respectively. The solution of the above equation will give us the 
answer we are seeking as follows: 


4 7T<p(x) = - I 

Jr 

-L 

L 


<Pn(y) 


8 (/) dy - V 


f 


<p(y)n 


8 (/) dy = 


r , , n 

8 (/) dy - (p(y)n- V Y 

r Jr 3 


V r 


1 5 (/) dy = 


(73) 


ip n (y) ipcosQ) 
+ r — \dS 
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We have used the result of Example 7, i.e., eq. (13), twice in the last step. The angle 
6 is between the local outward normal to the surface and r - x - y, see Figure 8. 
This is the well-known Green’s theorem for the Laplace equation. Note that a bonus 
of the imbedding method is that we know that by construction this equation will 
give <p = 0 inside the surface, a fact not at all obvious in classical derivations and 
that must be established separately. 

2- The Classical Kirchhoff Formula For the Wave Equation (The Exterior 
Problem)- We will now give the derivation of this well-known result based on eq. 
(71), i.e. using the imbedding method. Referring to Figure 9, we imbed the prob- 
lem in R 3 and use the fact that M n = 0 in eq. (71) to get: 

U 2 f = -<p n 6(f) - v-[<pnd(f)] (74) 


X 



(Kirchhoff) data surface 

Figure 9- Problem description for deriving the classical Kirchhoff formula, 
/ > 0 outside the surface and n - v / on the surface 


The Green’s function of the wave equation in the unbounded three dimensional 
space is: 


G(x, t, y, t) 


0 T > t 

6(T-t + r/c)/4-nr t t 


(75) 


where r = \ x - y \ . Here (x, t ) and (y, r) are the observer and the source space- 
time variables, respectively. We usually use the following symbol: g - t - t + r / c. 

We are first interested in the solution of the equation: 
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□ 2 = q (x, t) 6(f) 

The solution of eq. (76) is 8,9 ’ 29 : 


4 n if/( x, t) = 


f 9 *’ 


6(g) 

t ) 6(f) —— dy dr - 


L 


q(y,t-r/c) 


r jf = o r 

The limits of the integral in this equation are given below: 



Therefore, the solution of eq. (74) is: 

ip n (y,t-r/c) 


4 n <p (x, t) = - I 

Jf = 0 

-I, 


V, 


7 = 0 r 
W n (y,t-rlc) 


L 


n<p(y,t - r/c) 


>f = o 


f " 

Jf = o 


V, 


7 = 0 r 

( (f (y, t - r / c)\ 


(76) 

(77) 


(78) 


(79) 


Taking the partial derivatives with respect to the variable x of the integrand of the 
second integral, we obtain the following result: 


4 7l(f)(x, t)= I 

jf = 0 


C 1 cos Q (p t (y ,t — r/c) — cp n (y,t — r/c) 


+ 


/ ■ 

J/ = o 


cos 6 if (y, t - r / c) 


(80) 


dS 


The angle 6 is between the local outward normal to the surface and r - x - y, see 
Figure 9. This is the well-known classical Kirchhoff formula 34-36 . Note that we 
know by construction that - 0 inside the surface. This fact is not obvious in the 
classical derivation and, as in the case of the classical derivation of the Green’s 
theorem for the Laplace equation, must be established separately. Notice also that 
in our derivation we did not have to resort to the four-dimensional Green’s theorem 
for the wave equation. This is a great advantage of the imbedding method over the 
classical derivation. 


3- Radiation into Half Space- The Rayleigh Integral- We will derive two well- 
known results for radiation into half space by the imbedding method. Let the half 
space be defined by X 3 > 0. We are, therefore, interested in solving the following 
problem: 

□ 2 <^ = Q(x, t) 


(x 3 > 0 ) 


( 81 ) 
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<P (x, t ) = { 


(82) 


We will say nothing about the BC now because two possible BCs emerge from the 

imbedding method. Let us imbed this problem in R 3 in two ways as follows: 

ip (x, t) x 3 > 0 

-ip{x, t ) x 3 < 0 ’ 

~ . (Q(x,t) x 3 > 0 

Q(x,t) = \ 

1- Q(x,t ) X 3 <0 

and 

<P (x, t ) = { 


(Imbedding 1) 


ip (x, t ) 

ip (x, t ) 


^ . (Q(x,t) 

{ Q (x, t) 


(Imbedding 2) 


x 3 > 0 
x 3 < 0 ’ 
x 3 > 0 
x 3 < 0 

Imbedding 1- We can show easily that on the plane x 3 = 0 the jump 
A <p = 2ip{x\ t X2, 0+, t) = 2 q (x\^X 2 , t ), and also A (<9<£/<9x 3 ) = 0. Therefore, the 
extended function satisfies the wave equation: 

□ 2 ip - Q(x,t)-2q (x h x 2 , t ) d'(x 3 ) (x e IL) 

Using the free space Green’s function of the wave equation, we get 

Q(y,t-r/c ) C 26'(y 3 )q(y u y 2 ,t-r/c) 


(83) 


(84) 


4 7T<p(x,t)= f 

Jr 3 r 

Q(y,t-r/c) 


L 


+ 


X. 


r 

Jr 3 r 

d v2q(yi,y 2 ,t-r/c )n 


(85) 


dy 3 L 


dy 2 


J3=0 


When we take the derivative in the integrand of the second integral, we obtain the 
following result: 


Q(y,t-r/c ) 


4 7 T<fi(X,t)= f 

Jr 3 r 


+ 


f. 

Jr- 


( c 1 cos6q t {y\,y 2 , t - r/c) cos 6q (y\,y 2 ,t - r/c) 

+ 


( 86 ) 


dy 2 


Here we have used q t - dq / dt and the angle 6 is between the radiation direction 

x - y and the x 3 - axis. Note that the function Q is defined in eq. (82). This result is 
the solution of the following problem in the half space x 3 > 0: 
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U~p = Q(x, t ), (p{x i,x 2 ,0, t) = q(x i,x 2 , 0 


(x 3 > 0) 


(87) 


Imbedding 2- We can show that on the plane x 3 = 0, the jump A p - 0, and also 
A (dtp / <9x 3 ) = 2dtp(x\ j X2, 0+, t)/dx 3 = 2 g (xi 5 x 2 , ?). Therefore, the extended func- 
tion satisfies the wave equation: 

U 2 tp = Q(x, t) -q(x u x 2 , t)6(x 3 ) (xeR 3 ) (88) 

Using the free space Green’s function of the wave equation, we get 
Q(y,t-r/c ) r 26(y 3 )q(y h y 2 ,t-r/c) 


4 71 tp (x, t) = f 

Jr 3 r 

Q(y,t-r/c ) 


f 

Jr 3 


U 


r 

Jr 3 r 

2<i(yi,y2,t-r/c)-i 


(89) 


J r3=o 


<7^2 


Since the source in the last integral lies always on the y 3 = 0 plane, we can drop this 
subscript in integrand and we obtain: 


4 7T tp (X, t) = 


f 

J R ; 


Q(y,t-r/c ) 


f 

Jr- 


2q(yij 2 ,t-r/c ) 


< 7^2 


(90) 


Note that the function <2 is defined in eq. (83). This is the solution of the following 
problem in the half space x 3 > 0: 


□ 2 <^ = Q(x, 0, 


dp 

—— (xi ; x 2 , 0 , r) = g(xi,x 2 ,r) 

dx 3 


(x 3 > 0) 


(91) 


Equation (90) without the volume term is known as the Rayleigh integral 36 . 

4- The Imbedding Problem for a Second Order ODE- We consider the second 
order ODE lu(x) - u"(x) = k(x) for x e [a, b] c [0, 1] with two linear and 
homogeneous BCs which we denote as BC'i[m] = 0 and BC' 2 [u] = 0. Assuming 
that the solution exists, we want to find the solution by the imbedding method using 
the Green’s function of the Example 16. The situation is shown in Figure 10. 
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0 a b 1 


Figure 10- Imbedding an ODE problem for the interval [a, b] in the larger 
interval [0, 1] for the same differential equation 


Note that the ODE of Example 16 is the same as the one on the smaller interval 
This is required in an imbedding method. Let us extend the unknown function 
u{x) to the larger interval as follows 


ii(x) = 


( u(x) 
*0 


xs [a, b] 

, and 

xs [0, 1]\ [a, b] 


~ i k(x) x s [a, b] 

k(x) = l 

((3 xs [0, 1]\ [a, b] 


(92) 


See Figure 10 above. Notice that u satisfies the BCs of Example 16 given by eq. 
(62). This is required because if this condition is not satisfied, then we cannot use 
the Green’s function of Example 16 given by eq. (63). Note that we have 

lu = u" = f. We next find the following result valid on [0, 1]: 


lu(x) - u\x) = 

u"{x) + u{a ) 5 (x - a) - u(b) 5 (x- b) + u (a) 5'{x - a) - u(b ) 8'{x - b) = (93) 

k(x) + u\a) 8{x- a) - u\b) 8{x-b) + u (a) 8\x - a) - u(b) 8\x - b) 

Now let g(x, y ) be the Green’s function of Example 16 given by eq. (63), then we 
can find a representation of the solution as follows: 

ii(x) = I k(y) g(x, y) dy + u\a) g(x, a) - u\b) g(x, b) + 

Jo 


dg(x,a) dg(x,b) 

u ( a ) u{b) 


dy 


dy 


k(y) g(x, y) dy + 


, , dg(x, a) dg(x, b ) 

u ( a ) g(x, a) - u ( b ) g(x, b) + u (a) — u{b) 


dy 


dy 


( 94 ) 


FF&MKM- GFs Pt 1 (11-09) V5.nb 


41 


We are not finished yet. We can show that using eq. (94), u{a-) - 0 and 
u(a + ) = u(a ) give identical relations among the boundary values (BVs) of the imbed- 
ded problem. Similarly, we can show that u(b+) = 0 and u(b-) - u(b ) also give us 
another identical relations among the BVs of the imbedded problem. We also have 
two relations supplied by the BCs of the imbedded problem. These four relations 
will give us the four unknowns u(a), u(b), u'(a), and u'(b) which we can then substi- 
tute in eq. (94) to completely solve this problem. 

Looking at this procedure we note that all we have required of the behavior of the 
extended function ii(x) is that it should satisfy the BCs of the problem whose 
Green’s function is available in the larger domain. This means that any solution v(x) 
of Example 16 can be assumed in the extended region. Assuming that we have 
/ v(x) = ,/U) , then eq. (92) must be modified as follows: 


u(x) x e [a, b] ~ t k(x) xe[a,b] 

, and k(x) = l (95) 

v(x) xs [0, 1] \[a,b] 1 j(x) xs [0, 1] \[a,b] 

and eq. (93) becomes 

lu{x) = u\x) = 

u"{x) + [ u\a ) - v'(a )] S (x - a) - [ u'{b ) - v'(b )] 5 {x - b) + [u (a) - v(x)] 

6'(x -a) -\u ( b ) - v{b )] 8'{x - b) = (96) 

k(x) + [, u\a ) - v'{a)\ 8{x- a) - [ u{b ) - v'(b)] 8(x-b) + [ 
u (a) - v(x)] S'(x - a) - [u ( b ) - v(b)] 8'{x - b) 



Again using the Green’s function of Example 16, we find the representation of the 
unknown function as: 


f 


u(x) = I k(y) g(x, y) dy + [u{a] - v'(a)\ g(x, a ) - [u{b) - V{b)\ 

dg(x, a) dg(x, b) 

g(x, b ) + [u ( a ) - v(x)] — [u (b) - v(b)] 


(97) 


dy 


dy 


Two relations are obtained from u{a+) = u(a) and u{b-) = u(b ) among the BVs of 
the imbedded problem. Together with the two BCs of the imbedded problem, we 
can then find u{a), u(b), u'(a), and u'(b ) and substitute in eq. (97) to find the solu- 
tion of the imbedded problem. 
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The implication of this result for wave propagation problems can be substantial. It 
means that we do not have to assume that fluid inside the moving body is at rest 
with the conditions of the undisturbed medium in the derivation of the FW-H equa- 
tion. We can take any solution of the wave equation U 2 <p - q inside the surface and 
then apply the imbedding technique. This will allow us to use the free-space 
Green’s function and perhaps not have some of the shortcomings of the FW-H equa- 
tion such as the cancellation between the contributions of the surface and volume 
integrals. 


4- Concluding Remarks 

In this paper we have discussed some applications of GF theory principally to wave 
propagation problems. We have first given a fairly complete tutorial on GF theory 
and then have utilized most of the properties of GFs in the examples and applica- 
tions presented. We have emphasized generalized differentiation here which is a 
very important operation on GFs with many applications. These are some of our 
main conclusions: 

1- GF theory is a very powerful analytic tool in physics and acoustics with a wide 
range of applications, 

2- The derivation of many classical results in physics and acoustics can be enor- 
mously simplified by using GF theory, 

3- The power of the Green’s function technique is substantially increased if GF 
theory is utilized as in the imbedding technique, and 

4- It is possible to use functions other than the null function inside a moving surface 
in the imbedding problem used for obtaining the FW-H equation. Some of these 
choices may have advantages over the use of a null function in the extended region 
while still allowing one to use the free-space Green’s function. 

Generalized function theory is a very powerful tool in mathematics that every physi- 
cist and mathematically minded engineer should strive to learn. We hope that we 
have managed to impart to the readers our enthusiasm in learning and using GFs in 
their work. 


5- Notes 



FF&MKM- GFs Pt 1 (11-09) V5.nb 


43 


Note 1- To understand GF theory deeply, one must have some knowledge of 
Lebesgue measure and integration theory. Unfortunately, this theory is abstract and 
requires some mathematical maturity to master. Riemann integration is taught in 
undergraduate calculus courses and the students of physics and engineering see 
many applications of Riemann integrals in their work. On the other hand, one rarely 
has to evaluate the Lebesgue integral of a function but the properties of Lebesgue 
integrable functions are essential in solving many problems of physics and engineer- 
ing. It is important to remember that Lebesgue integrable functions constitute the 
largest and wildest class of functions ever encountered in applications. This prop- 
erty alone makes this class most suitable for use in applications. Lebesgue integra- 
tion theory gives us the definition of the integral, the properties and the theorems 
needed to work with this class of functions. Furthermore, this theory sheds light on 
differentiability and continuity properties of functions in the most general setting. 

To learn about Lebesgue integration, one starts by defining the measure of a set. 
For sets on the real line, the measure is basically a generalization of the notion of 
the length of an interval. Sets for which we can define a measure are called 
measurable sets. One can safely assume that he/she will never encounter a set that is 
not measurable in any application because defining a set that is not measurable is 
very difficult and involves the use of esoteric nonconstructive techniques or the use 
of nonstandard analysis. Next one defines a measurable function on the interval [a, 
b 1 as a function fix) for which the set {x \ fix) > c, c e R} is measurable. Finally, 
the Lebesgue integral of a measurable function is defined as a significant 
generalization of the Riemann integral leading to the class of Lebesgue integrable 
functions L which is a proper subset of the class of measurable functions. A 
function / is Lebesgue integrable if and only if | / | is. The definition of Lebesgue 
integral is quite abstract and subtle but one can become very comfortable with the 
use of the powerful results that are subsequently obtained. The extension of the 
Lebesgue integration theory to multidimensional space is feasible and fairly 
straightforward. 
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It is important to know that all finite and countable sets are of measure zero. There 
are uncountable sets that are of measure zero, e.g., the Cantor ternary set on [0, 1] 
obtained by dividing this interval into three parts and removing the middle third 
open interval and repeating this process ad infinitum on the remaining closed inter- 
vals. If a property holds on a set S with the exception of a set with measure zero in 
that set, we say that the property holds on S almost everywhere. For example, we 
can prove that two measurable functions that are equal almost everywhere on the 
same interval have the same Lebesgue integrals. That is why there is no regular 
function that behaves like the Dirac delta function because if there were, its 
Lebesgue integral would be zero (remember that the Dirac delta function is zero 
almost everywhere and, therefore, viewed as an ordinary function, its Lebesgue 
integral must be zero). A very important theorem that one should know is that a 
bounded function fix) on the closed and bounded interval [a, b ] is Riemann inte- 
grable if and only if it is continuous almost everywhere. 

One should realize that currently several different approaches exist for learning 
Lebesgue measure and integration theory. We are not familiar with a particularly 
superior approach. We believe that one requires a concentrated and dedicated effort 
to master the subject using any approach. There are many excellent books available 
on the subject. Two well known and lucidly written books are by Riesz and 
Nagy 37 and by Hartman and Mikusinski 38 . A recent excellent textbook is by David 

M. Bressoud 39 . This book covers all the necessary background as well as advanced 
subjects on Lebesgue measure and integration that one needs in applications. See 
the references in this book to many other books on this subject. 

We conclude this note by giving here Little wood's three principles that sum up 
nicely several important conclusions of Lebesgue theory 39,40 . In Bressoud's words, 
these are: 

i- Every measurable set is almost a finite union of open intervals, 

ii- Every measurable function is almost a continuous function, 

iii- Every convergent sequence of measurable function is almost uniformly 
convergent. 
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The word "almost" here means that the statement is true except on a set with arbitrar- 
ily small measure. Bressoud describes how these principles can be used to extend 
theorems proved using open sets, continuous functions and uniformly convergent 
sequences to much more general situations (Bressoud 39 , page 192). The second 

principle is actually Lusin's theorem and the third is basically Egorov's theorem 39 . 
Nikolai Lusin, who was a student of Dimitri Egorov, was the founder of the 
Moscow School of Mathematics known also as Lusitania. For those interested in 
the history of mathematics, a recent book by Graham and Cantor 41 gives the (often 
tragic) history of many brilliant Russian mathematicians of the twentieth century 
involved in the development of modern mathematics with amazing depth and 
breadth. 

Note 2- To the best of our knowledge this notation was introduced by F. Farassat in 
mid-nineteen seventies 42 . It was routinely used for many years in the course on 
partial differential equations ApSc-215 of The George Washington University by 
this author. Ram P. Kanwal gives credit to Farassat for the introduction of this nota- 
tion into GF theory in the first edition of his book on GFs 11 . This notation seems to 
appear most often in the literature of aeroacoustics. 

Note 3- Equations (37) to (39) were first published in these compact forms by F. 
Farassat in 1977 43 . It must be mentioned that several variations of these results are 
available in some books on partial differential equations but not in the simple forms 
presented here. 
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